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ABSTRACT 



A brief eunmry of the literature pertaining to the theory of 
cueuea la preaented to acquaint the reader vdth the mathematical 
techniques available. 

A quantitative analysis of the landing process around an air- 
craft carrier is made with a view towards a more efficient utiliaa- 
tion of the landing platform. The theory of queues supplies the 
mathematical techniques, thus making it possible to come to gripe 
analytically with the problem of reducing the landing time. 

The problem has been idealieed severely but 1s applicable under 
certain conditions to the actual landing process. The results may 
be used for investigating the distribution of the means of (a) the 
aircraft waiting time, (b) the time wasted by the aircraft carrier 
during the landing operations, and (c) the number of aircraft waiting 
provided we know (or assume we know): 

(i) the minimum safe landing interval, 

(li) the distribution of the landing times (landing rate), 

( 111 ) the distribution of the time intervals between aircraft 
entering the landing pattern (departure rate). 

If we l^ve not assumed the above, this information must be based on 
data obtained from, observation of actual landing operations. 

A method is also indicated whereby the problem may be analysed 
by a quasi-en^erical technique. This is a >veU known substitute for 
an analytical treat^nt if the mathematics become intractable. 
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V/ith the advent of the jet age, terminal congestion eind resulting 
aircraft delays are of Increasing ln 5 >ortance In 1 branches of avia- 
tion, The scientific approach provides a means for predicting how 
these delays yrlll be affected by varying condltiOTis, and so makes possi- 
ble an economic valviatlon of efforts directed towards reducing these 
delays. 

Operational research methods have a threefold application in this 
connection; (a) in ascertaining the wide field of practical problems 
susceptible to treatment; (b) in obtaining the basic data needed to 
enable the delay problems to be specified; (c) in selecting approp- 
riate solutions, 

oueueing theory supplies certain mathematical techniques ■vdxereby 
the operations analyst can come to grips analytically Tilth the problem 
of aircraft delays. Even where mathematics become intractable, there 
are good prospects of empirical solutions by means of electronic 
calciiLators . 

The purpose of this paper is to apply queueing theory to a prob- 
lem of seme inportance In naval aviation, carrier aircraft landing de- 
lays, It folloTra that the results may allow a more efficient utill- 
aation of the landing platform. Two things are of particular interest 
in the general queueing problem, the else of the q\»ue (stack) and the 
Tiaitlng time of the customer (aircraft delay). 

Chapter 1 gives an analysis of the general queueing problem and 

lU 



certain technical implications vjhich, while they ar« not strictly 
integral to the problem, throw considerable light on the report. Chap- 
ter 11 presents a detailed analysis of one aspect of the carrier air- 
craft problem, a siimmary of the results, the conclusions, and further 
considerations of statistical implications. 

This study was undertaken at the Iftiited States Naval T ostgradviate 
School during the latter half of the academic year 1953“54. -’he writer 
is deeply indebted to Professor Charles C. Torrance whoso constant 
criticism and enco\irageiDont were invaluable. To him go my foremost 
thanks. The tvriter is also indebted to Professor A, Boyd Ifewbom for 
reading the paper and for naking valuable suggestions and remarks. 
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QUEUEING THEORY 



IJL Introdxictlon to the Queueing Problem. 

The study of queues has been of interest to matheraaticiajis and 
engineers for the past forty years. <^ueueing occurs vdien demands for 
service are greater than availability of service, and It becouBS neces- 
sary to postpone the demands by a system of queueing or marshalling. 

’Ixile queueing is a natural process, marshalling is an atten^jt to do 
something about it. f^ueueing, however, is essentially a tenqsorary 
phenomenon, as othervd.se the queue would grow indefinitely in size. 

Thus one of the first problems in any queueing situation is to determine 
whether or not the demand will outstrip the service mechanism. V.hen 
this does not happen, so that the qxieue is continually returning to 
zero size, then the situation Is said to be in equlllbrivun, or, in statis- 
tical terminology, the process is stationary. Vhen the demand equals 
the service capabilities the situation becomes very delicately balanced; 
a slight x^uction in demand yields a sensible queueing system in eq\ii- 
librlxxm, while a alight increase results in an ever lengthening queue. 

The theory of queues has a special appeal for the scientist interes- 
ted in stochastic processes, for xmder certain conditions it provides 
an example which is both stationary and non-I/nrkovian. Thus one may 
expect in studying the queue to gain insight into other stochastic 
pAienomena, and to acquire a valuable facility in the handling of the 
relevant techniques. In addition, the theory has an astonishing range 
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of applicationo. A ctirrent liraitatlon of th« theory is that the nun^- 
ber of papers that deal vrith the theoretical side of the subject la 
small. And unfortunately those that exist are rather heavy reading, 

1.2 '^ueueing system. 

The s impl e nueuelng system is concerned not so imch vd.th the 
matching of supply and demand, but with a series of customers demand- 
ing seivice at a sin^e counter and waiting in turn to be served. Cus- 
tomer is here used in a technical sense; in the present paper it should 
be ecuated with aircraft. To complete the specification (Kendall 2 ) one 
must give a careful account of: 

(a) the input -process, 

(b) the qxieue-discipline , 

(c) the service-mschanism, 

1,2,1 Input-process, 

The input la random whenever the chance or probability that a 
customer vdll arrive at any Instant of time is the same for all instants 
past, present, and future, and also is unaffected by any past or future 
arrivals. The random input is often called a "poisson input” because 
of the use of the PoissoQ distribution to describe the arrivals, A general 
independent input (Llndley 3 ) is one where the intervals between the 
arrivals of successive customers are independent with identical pro- 
bability distributions, but this distribution is arbitrary. The input 
is regular when the cxistorasrs arrive at a constant rate. The input is 
correlated whenever the intervals between arrivals are not Independent, 

1.2.2 ^ueue-<iiscipllne , 

The cueue-discipllne is the rule or moral code determining the 
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nisnndr in which the customers form up into a queue, and the manner in 
vAiich they behave v;hile itaiting. In the simplest case they line up be- 
fore the single counter and airait their respective turns. The theory 
of multiple queues or many servers seems, etxcept under simplifying as- 
sun 5 )tions vdiich do not alv^ays correspond to reality, to be a problem of 
considerable difficvilty and x’dll not be considered here. 

1,2*3 Service-mechardsm, 

The service— mechanism may be described as the mode of service. In 
this connection the service time is of basic Importance and is the time 
that elapses while a particular customer is being served, 

1.3 Traffic density, 

k parameter of special importance is the traffic intensity or den- 
sity € • defined as the ratio of the average arrival rate to the average 
service rate. This may also be expressed as the ratio of the mean ser- 
vice time to the mean arrival interval, 

1.4 Itarkovian properties. 

Before the mthemtical theory is presented, an interesting quali- 
tative feature of the problem should be commented on, 

Suppose that the state of a stochastic system at time t is described 
by a rsindom function X(t)j the stochastic process is then said to be of 
Markov type (Doob 9) if a. knovtledge of the present value of X(t) makes 
all information about its past history irrelevant to a prediction of its 
future behavioior. Such a process may become non-Iiarkovian if part of 
the information contained in X(t) Is suppressed, for a knoid.edge of the 
previous behaviour may enable some of the suppressed information about 
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the present state to be reccTesTed* The queueing process with a Pois- 
sonian input and a single server is Ifeirkovlan if the present state of 
the process is described by the pair of random variables q and v^, vdiere 
q is the Instantaneous queue size and v^ is the expended service time 
of the customer at the head of the queue) in general it ceases to be 
f.Iarkovlan if the state of the process is measured by the queue size 
alone (Kendall 2 ). The only exception to this statement occurs v»hen 
the service time has a negative exponential diatidbution, 

1.5 Poisson input. 

Kendall ( 2 ) gives a general review of some points in congestion 
theory. A presentation is made of the PoUaczek "equilibrium” theory 
for the single counter queue fed by an input of the Poisson typo and 
associated ivith a general service time distribxition. It is pointed out 
that although the stochastic processes describing the fluctuations in 
queue size is not in general Markovian, it is possible to work instead 
with an enumerable liarkov chain if attention is directed to the epochs 
at Tdiich individual exist omers depart. The ergodic properties of the 
chain are Investigated with the aid of Feller’s ( 8 ) theory of re- 
current events. 

The simplest hypothesis about the input is one xvhich states that 
the customers arrive "at random", the number of arrivals in time T 
being a Poisson variable of expectation T/a, where a is the mean 
airival interval. The time interval t between two conseexitive arrivals 
will then have the negative exponential distribution 

d F (t) = e ^ , o < t < oo . 

o. / 

This is seen if one recalls that dF(t) is simply the probability of 
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no arrivals in a tins interval t. The successive t-variables vdll be 
statistically independent. This is the hypothesis adopted by Kendall, 
1,5*1 r«an waiting time vri.th general service time distribution. 

Let the queue be of the type vdth a single server fed by a I ois- 
sonian input and let the traffic intensity 6 be loss than unity, so 
that the system is not saturated. Let q be the size of the queue which 
the departing customer leaves behind him, not including himself, but 
including the next person to be served (n may be zero), '^’he next per- 
son to bo served has a service time v [t(v) - bj , say, and suppose 
that during this time r new customers anrive, "'hen conditionally r is 
a r oisson variable of mean value v/a[ and v has the service time dis- 
tribution, Lot q* be the size of the nvieue vdiich the next person leaves 
behind him. 

Assume that statistical equilibrium exists. Then the random var- 
iables q and q* have the same marginal distribution, and in particular 
their mean and mean sqixare values must be equal and finite; 

'• rcVi < ) e . 1,2 

The variables q and q' are related by the fundamental formula 

1,3 

or in terms of the S notation, 

%'= 'll-' 1./. 

v/here ^ is zero for all non-zero q, and ^ (O) ; 1 . It is 
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iraportemt to note the following conaeriuencea of the definition of the 
function S ( %): 



/ dnci (c^X I - O - V- ^‘5 

On forming the expectations of both sides of Eq, 1.4, it will be 
found that 

E(S)' i-e(-) •- I - ‘=4 ■ I- € . 1,6 

This is the chance that 5 is non-zero; that is, it is the probability 
that an incoming customer will not have to mit. On squaring both sides 
of Eq, 1,4 aixi making use of Sq, 1.5, one finds that 

— -t A-'-t 

= ^ i (Z A.'- I ) . 1.7 

Take expectations and note that r is independent of q and ^ , It 
foULows that 

£(c^) E (/^) + <T [n.i>tr -03 

2 ^ I - £( M i 

s ^ VCLn>. V 

2<xC<X-b) 

Nov/ suppose that a departing customer leaves q custcraers behind 
him, and let his own waiting time and service time be respectively w and 
V. Then q is the number of arrivals in a total time w + v. Thus, 

e (V) = . 1., 

o. 

Using Sqs. 1.8 and 1,9 one finds that 
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This is essentially Kendall's resxilt. It is convenient to express 
the resvilt in the form of the two ratios, for in applying the results to 
aircraft delays, one usually normalizes time in units of the mean arrival 
inteival or mean service tins (minimum safe landing interval), 

”/hen the mean airrival interval and the meain service time are equal 
to each other (€■!), the expected waiting time is infinite, thus show- 
ing that the intuitive solution € » 1 is not the best result. 

From Eos. 1,10 or IJLL one can see that the minimum expected wait- 
ing time will be obtained if the service time distribution is of the 
form 

Ftv) . O C''<k); FCv^=l Cviw) ^ 1.12 

that is, the service time is constant. On the other hand, if the service 
time is of the form 

. ^ . 

dF(v)= e dh/. . o < V < «o 1^3 

b 

the expected waiting time is twice the minimum value » l). 

It is tempting to think that the expressions for E(q) and E(w) may 
be valid for more general input processes; but this seems to be a fallacy. 
The difficilty is that r and q are no longer statistically independent. 
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Kendall shows that the Markov chain associated vrLth the queueing 
process is irreducible and the states are aperiodic. A further classifi- 
cation of the states depends on the value of the traffic intensity, € • 
V.hen € is less than unity, all states are ergodic; YJhen € ia greater 
than unity, all states are transient; in the critical case when € is 
equal to unity, the states are all recurrent and null. (Feller 9 ) 

1 . 5,2 Aircraft delays, 

Bowen and Fearcey ( 6 ) give an elementary analysis of the effect 
of control procedures on the flow of air traffic, v/hlle Fearcey ( 7 ) 
presents a more rigoroxis treatment of the problem and derives the dis- 
tribution of delays under different traffic conditions. The hypothesis 
of random arrivals is used, so that one is led to the negative exponen- 
tial distribution for the airival intervals. See Eq, IJ., 

The assun^Ttion is made that the miiiiniun safe landing interval 
(service tine) is constant. Thxis we have Kendall's solution with a con- 
stant service time in a different notation. 

Fearcey '3 treatment is based on a rather congxLlcated procedure of 
successive integrations but once the assun^Ttion is made that the system 
is in a state of statistical equilibrium the analysis must lead to the 
same results, Cne feature of interest is the derivation of the proba- 
bility of successive delays. By successive Integrations Fearcey finds 
that in general the probability of n-1 aircraft being successively 
delayed is (stack of n-1 aircraft) 
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P(y>) 

M'r.) I - ^ - n.)^ /\tr\-rx\ 

"TTi — ' 

Mo] z 1 ^ 

T-. 

Bell ( 5 ) gives a someivhat more general treatment by using a moment 
generating procedure. 

1.6 General independent input. 

The center of interest of Lindley's ( 3 ) PQP*r is the ivaiting times 
of the customers. A relation betvreen the miting times of successive 
customers is found, and this makes it possible to calculate the wait- 
ing time distribution of any customer. 

This relation is also used to obtain necessary and sufficient con- 
ditions for the existence of a stationary (equilibrium) state. The paper 
concliides Avith an investigation into the stationary state in some spec- 
ial cases; particularly in the case Avhen the customers arrive at regriLar 
intervals , 

1.6J- Jaiting time distribution. 

The ''ueue here is again of the type vd.th a single server but vdth 
general service time and input dlstribvttions ; the interirals betA'/een 
successive arrivals may have any suitable distribution. 

Let tj. be the time interval betvreen the arrival of the r and 
r + 1 customer, and let s^ be the service time of the r customer. 

The following assAimptions are made: 
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1. The are independent random variables \vith identical 
probability distributions and the mean, E(t ), is finite; 

oO 

E (t^) = j cLt^ < OD 

o 

N) 

2« The s^ are independent random variables ■(vith Identical 
probability distributions and the mean, E(s ), 1s finite. 

T 

Thus the two sets of random variables and [ (r • 1,2, ...) 

are statistically independent. 

Let be the waiting time of the r^^ customer. It is evident 
that the following relations exist (Fig. lU.): 



r 



t'. 




a 






/ aaiMval ) 
^ -t 1 rn t / 



^ ^ \ *1 " < ^r ^ »r ) 

1 (tj. i ^ ♦ Sj.j see tjJ.) 



r + 1 






1.15 



Let ^ ■ ®j. ■* 'V independent random variables 

(by virtue of the preceding assumptions) and have identical probability 
distribxitions ; E( | u^ | ) is finite and u^. is Independent of wij,. Then 



- IV ♦ Ur 


(Wj. + Uj. > 0), 


1.16 


- 0 


(Wj. + u^ 1 0). 


1J.7 
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Nonv if wo aissixnB that the first customer does not hetve to wait 
(wj^ • 0) we have the follovrlng relations: 



w - o • 
• ) 


w^-- u. • 


II 






W : 

r\^-i 


r\ rv- 1 


• • 


. , + u, 


1.18 



and the unexpected fact develops that the waiting times depend on only 
the difference between the service time and interval time. If we define 
F (x) as the probability that w_ i x, and G_(x) as the probability that 
u i X, we can sli^tly modify Lindley^s procedure and find the distri- 

T 

st 

bution of waiting times of the r + 1 customer by using Eq, 1.13, The 

u are independent statistically. Thus the distribvrtion function of 
T 

Fj. + 1 is the convolution of u^ ♦ u^^ + ..... + (Doob 9) and 



pK„ i X) = G. Oi.) X do,. (V.) 

Using the definition of Gp(x), to have 



1J.9 



: p - X ; X ; x). i^o 



Eq, 1 J.9 establishes a recursive relation such that the individ- 
ual waiting time distributions may be found. To find the final behavior 

of the Queue , we must find the limit of F .. (x) as the number of cus- 

r+1 

tomers ccntinually increase. Let E^, be the event: 



U t 



U. t K • 
) 






i X 



- X 
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Noiv the sequence of events ^Eyj decreases to a limit E as r approaches 
infinity, that is: 
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= L*n\ pCErv) : plE) 

ft —V o® ' 



1,22 



by a theorem of probability measure. In the limit, then, 



f FC y-u.) d G (u.) 



1.23 






= / FC^)dGCx-'^) 
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This integral equation is difficult to solve, but it is possible to do 
so once G(u) is specified and certain conditions are satisfied (Paley and 
Weiner 11-). The solution of either integral equation gives the dis- 
tribvttion function of the waiting time, that is, the probability of 
waiting a time not more than x. The distribution function for any pai*- 
ticular customer can be calculated by performing the necessary inte- 
grations; but this is also difficxilt to do. 

Lindley uses Feller *s theory of reoxrrent events and the strong 
law of large numbers to investigate the behavior of F(x). 

The event of "not having to wait”, called event A, is a recurrent 
_ event, since the qvieueing process may be considered to make a fresh start 
whenever a customer does not have to wait, and its previous history is 
quite irrelevant to subseqvent developments. 

To investigate the behavior of F(x) one must examine the relation 



E ^ e (S) - E ct ^ 



1,25 



as follows: 



E(u) > 0: The strong law shows that the event A is transient and 
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that with probability 1 there ;rUl be soxoe point in the quetieing process 
where every subsequent customer will have to wait, 

E(u) < 0: The strong law shows that event A is certain with a finite 
mean rectirrence time. It should be dear that it is not periodic. 

E(u) - 0; The event A la certain, but the mean recutrence time is 
infinite. Hence it is certain that there will always be another occasion 
on which a customer does not have to wait but the probability of any 
particular customer not having to imlt tends to zero. 

Thus a necessary and sufficient condition that the waiting time dis- 
tribution function tends to a non-degenerate limit, as the number of cus- 
tomers increases, is that E(u) < 0 or u - 0 certainly, 

Lindley sdves the integral equation when the arrivals are random 
that is, 



These distributions are the Type HI curves of Pearson conmonly called 
the gannia distributions (Mood 10), and are of the form 




1.26 



and also for 



d & (^) = X e ^ EC'^) : 



_2 

A 



1.27 




1.28 



with the given parameters n and X . 

And as expected, the probability of not having to wait is 



1 - e . 



EC-t) 



1.29 
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The equation to be eolved for regular arrivals, where the arrival 



interval is the unit on the time scale, is 



Fix) = d&, : 



1.30 



with the sejTvice time distribution as given by Eq« 1.26, 

The solution is difficult to develop, but many may be expressed in the 



form 



/Z ; X 

. F U) = < - C i e " 

L i I 



c. : C. 


3 nd 


t. \ 

' ■ 


1 . , 
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. n • 
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21- 

: - , , n -t- f 


:h . .. z. 


9 n d 
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) 1.32 

O (a = o, I,-- - V^) ^ 



The ustial form of service time distribution is one vAiose freqvjency 
function is zero at s increases to a maximum and then decreases, 
with a long tail, to zero as s approaches infinity. In some cases it 
may decrease steadily from some non-zero value at s ■ 0. To a good 
approximation these functions are represented by the gamma distribu- 
tions (n - 0, 1), 

The results of the calculations sho7f that regular arrivals are 
roughly twice as efficient as random arrivals, i.e., the probability of 
not having to wait for regular arrivals is roughly twice as groat as 
the random arrivals. Lindley's calculations for a few values of <£ are 
listed in Table 1 to show the nature of the results (xmits are arbitrary). 

U 
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Probability of not havirv; to wait 



ISsan vraltlng time 



n " 0 



€ 


Regular 

Arrivals 


Random 

Arrivals 


Regular 

Arrivals 


Random 

Arrivals 


0*830 


0.314 


0M7 


1.823 


4a66 


0.715 


0.511 


0.286 


0.684 


1.786 


0.635 


0.642 


0.375 


0.349 


1,042 


0.500 


0.797 


0.500 


0.128 


0.500 


n - 1 










0.830 


0.386 


0.167 


0.834 


3.124 


0.715 


O.608 


0.286 


0.285 


1.339 


0.635 


0.743 


0.375 


OJ.32 


0.781 


0.500 


0.884 


0.500 


0.039 


0.375 



Table 1 



1*7 Correlated Input. 

In view of the difficiities encountered v/ith the random and inde- 
pendent inputs, it should be readily apparent that the inclusion of con- 
ditional probabilities in the input process will make the problem even 
more difficult, Ae mi^t be expected then, this problem has not been 
solved analytically. 

The use of the Poisson input in air traffic loads to a complete 
analytical solution of the stack-delay problem, but unfortunately once 
aircraft are in the air, the input is not even approximately random, 
but has the form of a regular pattern with largo errors. It is virtually 
ijDi5)os3ible to find the actual distribution vdth any accuracy, as the 
amount of data required for an unbiassed san?le is prohibitive. The key 
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to the difficulty lies in the fact that aircraft in the air proceed on a 
predetermined schedule. This detailed information in the schedule should 
be regarded as a datum of the problem and used in its solution. 

The first significant contribution of the electronic calculator to 
the queueing problem has been its use in the study of the flow of 
scheduled air traffic. The Research Laboratory of Electronics at the 
Massachusetts Institute of Technology { k ) are the pioneers in this 
endeavor. They undertook a program to determine, as far as practicable, 
the quantitative relation between the degree of control of aircraft 
enroute to an airport, and the resulting congestion and delay. 

Actually there ai^ two possible problems involved in the control 
aspects of air traffic. One is the maintenance of a pre-assigned schedule, 
and the other is the reduction in the congestion around the terminal or 
landing strip. It should be noted that this is an unusual situation in 
that the aims of the custcaner and the server are the same: (a) to ensure 

safety in the airj (b) to allow the most efficient utilization of the 
landing platform) (c) to reduce the total time by avoiding long delays 
at the terminal. This is not always true in the general queueing problem. 
In the UTr report, IBM punched-card machines vrere employed to make 
a theoretical analysis of the resulting congestion at a single landing 
strip vdien aircraft are scheduled to arrive there in some proper se- 
quence, but fall to meet such schedules according to certain simple 
deviation statistics (rectangular, triangular, and p>arabolic). These 
deviations statistics are assumed to have a finite spread, that is, 
deviations in excess of a specific amount have zero probability. This 
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is within the bounds of reason from the pilot's viewpoint. Numerical 
results for the distribution of the resulting stack and total delays 
are presented. The former are con?)ared with those viiich arise from 
the more random Poisson arrival distribution. The amounts of conges- 
tion and stack delay are found to be considerably smaller under the 
new conditions, particularly rihen the traffic is heavy, A significant 
range of the parameters has been covered by the numerical analysis, 
supported by strictly analytical methods idierever possible. Thus there 
is established a quantitative relation between the time-keeping errors 
of the aircraft enroute to an airport and the resulting terminal con- 
gestion and delay. 

The density of the traffic flow to an airport is described by the 
usTial traffic parameter € , defined as the ratio of the average arrival 
rats at the airport (number of planes par \mlt time) to the maxLraum 
allowable acceptance rate. This is also the ratio of the landing in- 
terval (service time) to the arrival interval. 

The authors found that it was convenient to use time in discrete 
units rather than cm a continuous basis. The ndniraum safe landing 
interval, t^, is chosen as the unit of time, and forma the base for 
both the normalization and quantization of all other times appearing 
in the problem, 

A sample size of 1,000 planes was selected for use in the numerical 
analysis. In the programming procedure, each plane in the sample is 
represented by a cai*d, which is identified by a number J, with J greater 
than or equal to 1 and less than or equal to 1000. Thus J defines the 
scheduled arrival order. The scheduled arrival times tj are then as- 
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signed using a random number table. 

Between its scheduled take-off and actual arrival times the aircraft 
is subject to a delay rj, which la assigned (pwnched onto the J^^ card) 
according to the enroute distribution used. Thus plane J actually ar- 
rives at time pj - tj "*■ number is also entered on the J^^ 

card. It is possible now for more than one plane to arrive at the 
same tinx and for them to arrive in a se'^uenoe idaieh differs from that 
scheduled. This leads to the formation of a stack of aircraft, and 
consequent delays idiUe vraltlng in the stack. It is then assumed that 
planes will be handled at the terminal in order of actual arrival, with 
the provision that simultaneous arrivals will be considered as ordered 
according to their originally scheduled sequence. 

After the preliminary details were comfleted, i.e., the t and r 

J J 

tables were prepared, and the initial conditions were chosen, the IBIJ 
machines ^vere prograunmed in detail for the specific cases, i.e., for 
each deviation distribution, a spread S, and traffic parameter £ . 

The principal results which cane directly from the machines were: 

(a) Frequency dlstribiitlons for the stack delays. 

(b) Average stack delay for each irun. 

(c) Cumulative flnequency distributions for the stack delays. 

(d) Fre''uency distributions for the total timekeeping errors 
(stack delay plus enroute delay), 

(e) Average timekeeping error for each run, 

(fy Cumulative fre-mncy distributions for the total timekeeping 
error. 
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V/ithout loss of generality, it vras assumed that all variation* 
were delays (simple change of scale). 

The ^•3^^1ts form a considerable source of data pertinent to a 
traffic of properly scheduled aircraft. The data are given in detail 
in the form of curves of stack-delay probabilities, prcgressive-delay 
probabilities, and average delays, for the three types of enroute-<ie- 
viation distributions considered. 

The curves give additional information regarding the effects of 
variation of the parameters. In particxilar, only relatively small 
differences in the results are caused by rather moderate variations 
in the shape of the enrotite-deviation distributions. The relatively 
small changes in the stack-delay distribution caused by making moderate 
changes in the shape of the enroute-deviation for given S made it 
posslKLe to average together the IBM results for all three basic shapes. 
The average curves consequently give stack-delay probabilities fca* an 
average shape enroxrt e-deviation distribution, which is sort of cross 
between the rectangular, triangular, and parabolic shapes. 

The effects of changes in the traffic peraia»ter € are practically 
the same for all three types of distribution. For values of € close 
to unity, the average delay changes fairly rapidly x-d.th € , while cor- 
responding changes in € , when it has smaller values, cause little 
change in the average delay. The staclo-delay distribution is also de- 
pendent in a similar way upon the traffic parameter, %on € is close to 
unity, the most pi*obable stack delay is almost equal to its average valve, 
bvtt as 6 decreases, the most probable delay rapidly approaches zero 
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with a amall average delay. The remaining variable parameter, the spread 
S, also has a decreasing effect vhLth decreasing values of £ « 

1.8 Tabulation of cases covered. 

The principal oases covered may be tabxilated briefly as follows; 

(a) Poisson input with general service time distribution (Kendall 2 ), 
Derives expressions for the expected waiting time and the ex- 
pected queue sise. 

(b) General Independent Input i.e., general service time and input 
distributions including regular arrivals (Lindley 3 )* 

Derives the waiting time distribution of any customer and the 
distribution function of the waiting time in the stationary 
case. 

(c) Correlated input i.e,, scheduled air traffic (Adler and Pricker 
U ). Use electronic oalciilator (lEU) to find the frequency 
diet rib\it ions for the stack delays, oxumilativo frequency dis- 
tributions for the stack delays, and the average stack delay 
for each run. The information is presented in the form of 

I 

useable graphs. 

1.9 Implications 

An interesting Interpretation of the MIT results may be vlsuallxed 
in the foUovdng manner. Consider the problem of bombing when a largo 
number of planes imy be directed to strike several targets. The usual 
practice, to avoid saturating the handling capacity of the airfield, is 
to stagger the formations in such a manner that all of the aircraft will 
not return at the same time. 
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For exarapLo, considar a 200 plane strike made up of 20 groups of 
10 planes par group. Now ea hare no idaa of the anroute darlation dis- 
tribution, but lat us use our imagination and say that the average 
anroute-deviatlon distrilwtlon described in Section 1.7 is applicable. 

Assume the expected lll^t time to target and back to be 10 hours 
for each group. Let the ndnljituin safe landing Interval bo approodlmataly 
2 minutes par plane so that it takes on the average 20 sdnutss to land 
a group (t^ ■ 20 ndnutaa). 

It is desired that the 20 groups arrive home in such a manner that 
the planned maximum stack is to be 2 giroups or less (20 planes). Larger 
stacks are not feasible due to lack of control facilities, safe altitude 
separation, etc. 

For a 10 hour fli^t, errors of * 20 minutes would not be deemed 
too large j in the notation of the IHT repjort, this will correspond to a 
spread of 40 minutes. If we normaliee the spread in units of t^, S - 2. 

If we apply the res\ilts of the 1£ET paper with € ■ 1, that is, with 
the departure interval also 20 minutes and with S ■ 2, we find the pro- 
bability of stacks greater than or equal to 2 groups is appioximately 
0.95 , In other viords, it is almost certain that moro than 2 groups will 
overlap at the field. If we reduce £ to O.S, that is, make the dop>ar“ 
turn interval 25 minutes, the probability that stacks greater than or 
equal to 2 groups will form is reduced to 0.06 . Thus it may be seen 
that quexieing theory, as developed by electronic computers ixsed in con- 
junction with good experimental data, has striking pjossibillties. 

This eocaxnple is of course trivial in that it is such a restricted 
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problem. But erren so, the results are well worth further study. The 
Poisson solution when € Is 1 has been stressed enough; for exanqple, the 
stack approaches infinity, and the ^waiting time approaches infinity. 

Thus there woxild be no solution in that event. VJhen € is O.B, the 
Poisson solution gives for the probability of 2 or more groups in the 
stack a value of 0.25 approrLmately, In conq^aring the Poisson solu- 
tion with the IBII solution we see the two results differ by a factor 
of 4> a considerable discrepancy. The key to the difficulty is the 
use of schedxile information as a datum of the problem. 
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C.^ILR .\IRGROT LANDI^D DELAYS 

2J. Fomulatlon of the problem. 

The prime desideratum in landing aircraft aboard a carrier is 
simply to land the aircraft as quickly as possible while observing the 
necessary rules of safety. The time required to land (and laimch) the 
aircraft aboard the aircraft carrier is of basic importance due to the 
extreme vulnerability of the carrier during the landing operations. The 
intuitive solution of the problem of minimizing the landing time is to 
mlnimiae the landing intervals by feeding the aircraft into the landing 
pattern as qulcldy as possible subject to the i^straint of minimu m safe 
landing interval. tVhile this is sometimes a satisfactory solution, yet 
sometimes the delays involved become excessive. Thus in attempting to 
minimize the landing tii», one must investigate the behavior of the air- 
craft diu*ing the landing process. This fact makes necessary a quanti- 
tative analysis of the landing process, with a view towaixis a more 
efficient xitilization of the landing platform. Thus it seems natviral 
to apply the theory of qwues in the analyeis. Let us consider the 
main aspects of the situation under the three headings introduced by 
Kendall to characterize a queueing process. 

As a preliminary step, consider the pattern of events around an 
aircraft carrier. All of the aircraft waiting to land are placed in an 
orbit above the carrier at a safe altitude tmtil the signal to begin 
lawiing is given. This may be vlsvialized as a ready-made stack (queue). 
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'.4ien the si^ial to land is given, the aircraft descend to the landing 
pattern to cconence their approach. There is a constraint in that onl 7 
a maximum number of aircraft are allowed in the landing pattern. Every 
time an aircraft lands it is replaced in the landing p>attem by one of 
the waiting aircraft. Thus one can visualise the single server quevie- 
ing process. 

First, there is the input process. In our context this is essen- 
tially the process by which the aircraft depart the orbiting stack (queue) 
and feed into either the landing pattern directly or into another stack 
(queue) of aircraft waiting to enter the landing process. The landing 
pattern may also be visualised as a queue with its associated piro- 
hlem of delays due to wave offs, etc. This is a problem for further 
study and is not considered here. Thus we may visualise two queues 
with a single server. It should be noted that the orbiting stack (high) 
is called a queue in a strictly limited sense. This artifice makes the 
analysis more understandable. The low stack (queue) is the feature of 
interest in that it is caused by some p^erturbation in the landing process. 
V/e shall assume that all the aircraft leave the high stack at a con- 
stant rate (continuous flow), i.e., at a succession of eq\iaUy spaced 
times, the difference between one time and the next being called the 

4t 

departure interval. Thtis the dep>arturo distribution is of the form 
(3(v)-0(n.<v)-GCv^zl (R>v)^ 

where r is the departure rate (constant). Second, there is the queue 
discipline, shall require that the aircraft land in the order in 
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which th«7 arrive at the low etack if one exiets. Third, there ia the 
eervice mechaniam, which is given by the frequency distribution of the 
landing times (serving time). In a limited sense we shall consider the 
landing pattern as a part of the service mechaniam« Landing time is 
here taken as the time from the cut signal given by the landing signal 
officer to the time the aircraft is safely forward of the barriers. VJe 
\7ill now make the drastic simplifying assungAion that the landing time 
is constant, i.e., the landing distribution is 

FCt) = 0 (Us) ; F (S) : I (U 

where 1 is the landing rate (constant). This assun^jtion is made for 
two reasons: (a) observational data ar<e lacking, and (b) it is desired 

to make the analysis tractable. It seems rather natural to assume that 
the landing time distribution is one Those frequency fioiction is exit off 
at the minimum safe landing interval. Increases to a maximum near some 
average value of landing time, and then decreases to another arbitrary 
cut off point as the landing time Increases* Some consideration will 
be given later to certain statistical in^lcatioos pertaining to the 
gathering of observational data. A second drastic assumption to be made 
in connection with the landing time is that \inder normal conditions 
(when the normal landing rate is in effect) the aircraft are landed 
aboard without ai^ delay for a certain length of time, and then the 
normal landing rate is changed abruptly, and continues at a new reduced 
constant rate for a period of time T^. The reason for this ireduction 
in rate may bo poor handling technique aboard the carrier, inexperienced 
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pilots, etc. During this period of tijne the low stack begins to form. 

It is now easily seen why we consider two queues in series with a sin^e 
server. 

An inyx»rtant parameter of the problem is the nvuiber N of planes 
initially in the high stack. Uhder normal conditions (when the low stack 
is empty) the tina for landing N aircraft is 
Tj. - N/r. 

Thus the problem is to investigate the delays under the two ccsiditiona, 

T > T and T i T^, 
o r or* » 

Before proceeding xvith the analysis it shovild be noted that the 
fundamental difference between the general q\jeueing problem and the pro- 
blem at hand is that we are not interested in the limiting case, that 
is, in an infinitely long quexwlng process. The stack (queue) assoc- 
iated with the cazrier piroblem has a finite bound N as a datum of the 
problem, and this must be Inclnled in the analysis. However, as the 
departure interval approaches the landing time, we mL^t expect from 
the work of Llndley that the situation is someidiat critical for short 
queues as iveU as for long ones. This is in fact the case. 

The notation to be used is as follows: 

X (t): landing rate aboard aircraft carrier (planes per unit of 

time). 

r (t): departure rate from high stack (planes per unit of time), 

€ s ^ : normal traffic density parameter also equal to the ratio 

of the landing time to the departure interval. 

a : landing rate reduction pjarameter. 

T ; total time a is diminishing the landing rate. 
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: total time during vrfilch non-aero lov/ stack exists. 

Tj. : total time for carrier to con^iLete landings under normal 



N 



conditions . 

number of planes initially in high stack 



S(t): number of planes in lew stack (low stack size). 



ac 



total additional time required for carrier to con^iLote 
landings due to formation of low stack. 

•X : aircraft delay. 

D ; total aggregate delay (say in plane-cdnutes) equal to product 
of average delay (?) and the total number of planes delayed; 
also e^ual to the product of T, and S(t). 

Note that the traffic density € is not a function of time. 

The numerical quantities of primary interest wHl be: (1) the ad- 

ditional time required for the carrier to congxLete landings due to the 
reduction of the landing rate; (2) the average stack delay over the per- 
iod during viiich a non-zero low stack exists;, and (3) the maximum stack 
delay suffered by any plane. 

The calculations fall into t;vo major groups, defined by i and 



T i T_, The first case to be considered will be T. > T . Figs. 2J., 
2.2, 2.5, and 2,6 show the corresponding departure rate, landing rate, 
and stack height for prescribed valxies of a and € in the ranges 

0 < a < £ <1. 

\’le can asstune without any loss of generality that commences 
when t - 0. Prior to this .time, the landing process is normal, say 
( € s ~ i I ) and there is no low stack. It is not essential to 
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the analysis ^vhether the landing process Is in progress at t ■ 0, or 
Just commencing. 



Now at the beginning of tine the landing rate is reduced to aJi 
and, provided a <6 , a low stack begins to form at a rate (r - a X). . 
r-ai^ - j((€-a). 

It is essential to keep in mind that no stack forms at all unless a < £ , 
and this condition is therefore impLiolt throughout the derivations. At 
time T^, the departure rate from the hl^ stack becomes sero, and, in 
consequence the stack no longer increases bxzt begins to decrease at a 
rate a . After a further period of time (T^-T^), the aircraft carrier 
resumes its normal landing rate. As the hi^ stack has deajred, the low 
stack decreases vdth the normal landing rate ^ . 

2.2 Stack slee S(t) and stack time T.. 

At time the low stack hei^t is 

S (T,) :T. JK6 -«) ; 2J. 

S(T«) Is equal to S(T) . See Fig. 2J.. Now the stack begins to de- 

*’ niHX « 

crease at a rate , and at time T^, 

S (T„ ) t -ocl - tT. -T,. ) o( X 

: i ( t - cc Vt. ) . 

Now S(T ) is non-negative if and only if 
o 

£ > To . 
if this is not the case and 




the stack clears entirely prior to the end of tine T^. See Fig. 2.2. 
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Thus have two cases to Investigate: 



Case A; -i- - Xa. > I . 



Case S; lEs- ^ > I 

o<. 



Case A: 



A- > To >1 

= 7 ? 



The time taken to dear the stack left at tine T is 

o 



S (T.) 



n. ^ t - oC "7° \ 

T 



= -n. U - . 



2.3 



Now the time that a nonHBsro low stack exists Is 
Tj ^ -r» + Xv ( <i - ) , 

r\ 



Zf. ^ -i- > I 

o< 



Case B: To > € 

The stack height at time T^ is 

S Ct^ ^ A ( e - . 



2.4 



2.5 



In this case the stack begins to decrease at a rate a ; thus the time 
for the stack to dear is 

SCx ) , T., X C6 
“x X 04 H 

r V 2.6 

■- • 
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The time that a nan-xaro stack exists Is 

T,-- Clr'O 



2.7 



Average lonv stack S, 

Case A: ^ ^ I 

T“_ 



The average stack height by definition Is 

sir) _L S (t ) d t 

-Ts ^ 



. I 



J "t )! (£-=<) 






r\ 






and after integration and simplification 



S : i^'To ( »‘ <5)1 



2 [<: -*-T^ (l-»<)l 
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Case B: 



Xi. - i. > » 
■r. 



The average stack S is easily seen from Fig. 2.2 to be 

s - ^ 5 ( c- - »< ) . 



2.9 
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2,3 AT»rage delay 'T' « 

The deriratlon of the average delay by the usual process for com- 
puting averages is rather unwieldy. It is much more convenient to com- 
pute T from the ratio of the average stack height s(T^ to the average 
landing rate (V) or the average departure rate r(i), all arerages 
being taken orer the period during which the stack Is never aero. 
That T can very generally be oon 5 )uted correctly in this maimer may be 
demonstrated in the foil owing way* 



S(t) 



-t t+dt t, -Tj -t ^ 

Stack Height vs. Time Plot for General Case 

Fig. 2.3 

Consider a low stack coanencing at t - 0 and ending at t • T^ as 
shown in Fig, 2,3 « Assume the variation of S(t) as arbitrary but subject 
to the following restraints: 

sto) -- S O, 

foa 0<t.< ^ 

2JLL 
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In terms of the notation previously defined, it foUorrs from Eqe. 2.10 
and 2cU that 



s(o-- 

o 

O i -t i ~T^ 



ft (t ) "] ci "t 



2.12 



This of course puts certain general restrictions upon r(t) and (t), 
since they must lead to a stack subject to Eqs. 2.10 and 2.11. Chranting 
these general restrictions, Eq. 2J.2 is valid regardless of the detailed 
form of either r(t) or X (t). 

j ^(t idt 

o 

That is, for such a stack variation, the total number of planee depart- 
ing the high stack equals the total number landed. 

Consider the planes r(t)dt ;diich arrive in the time Interval be- 
tiveen t and t dt. They Trill have to rrait say, until t •* t^^ to land, 
because of the time inquired to dear the stack S(t)« It foUotra there- 
fore that t]^ is a function of t established by 

S(0 = J*'j2U) di 

“'t ’ 204 

Tvhich simply states that all the planes in the stack at time t must 
have landed at a rate (t) by the time t^. From Eqs. 2.12 and 2.14, 



Now since S(T^) • 0, 

r 

J a(t)dt . 



2.13 
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2.15 




\j\ ^ (t )*] 

j rv(,t ) cit 



cit = 




clt 




I 



nhlch defines as a single valued function of t. Th^ul as ml^t be 
expected, the delay of a plane Is a function of the time it arrives. 

Now by definition, the average delay ^ is the total aggregate 
delay D divided by the number of plana s delayed, id\ere 

D = ( t,-t) a(t ^ dt . ■ 



That is, the planes arriving in the tine interval dt are delayed a time 
(t^ - t) and the delays must be sunned over the period during which 
the stack is never zero. And of course the number of planes delayed is 
simply the total number arriving during the period T^. Thus 



'T 



J ' (-t ^ ci -t 



2.16 



or In view of Eq, 2.13» 






/, * ( t.* t ^ n. (_i: ) d± 



2J.7 



Now the numerator can be integrated by parts as foUorrs: 
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r^s 




j n.(t ) d t = 


■tj n-Ctldt 


o 


^ ... 



t r o 




By Eqa« 2.10 and 2 J.5, 

- 0 wh«n t - 0, 

t-j ■ T vrh*n t - T . 

i • 8 



2.18 



2J.9 



2.20 



This sluply states that the aircraft vrhlch arrive ivhen the stack is zero 
are not delayed at all. Using Eq. 2J.5 to eliminate 

j\(M dt 



from the last term of Eq» 2.18, ee have 



J 




O 



i\(i) ol t 





and by a similar process 

j t rxC't'ldt c T - j 

® o 0 ^ ® ' 



2.21 



2.22 



Thus 








2.23 



Now dt^ is merely functioning as a variable of integration in the last 
term of Eq. 2,23» so that 
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Using the rather obvious definitions 



S Ct ) = 



\ 



T-s 




SCtU t 




^ XM d't 







I f\ (.t \ d t 

O * 



2.25 



2.26 



2.27 



and Eqs, 2 JJ2 and 2 J.3, vw hare the desired alternative form of Eos, 
2J.6 and 2.17: 

n" - 






J 2 (t) 



; 



and 






^ Lt) 



S (^) 



n. ( t) 
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Case A: > "T o > | 

o< nr 
• ^ 

Using Eqss 2e4# 2e26, and obserring Fig* 2J., ws find that the ar- 
erage landing rate is 

X(t) : -r j dt + J Jidt 

° T, 

ej 

"T'o ( • " ^ 



2.29 
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Thus from Eqa. 2«2S, 2,29, and 2.8, the average delay is 



*- 



' o 
2 ^ 



[(!-^)U£ ■^'Z. ). eZi (i-t\ 
L T. -r^ ' 
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Case B: ~^o ^ ^ > 

"T cK 



From Fig. 2.2 and Eqe. 2.26, 2.28, and 2.9> the average delay la 



Z 
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2,4 Maximum delay . 

In order to find the maximum delay, the delay of planes arriving 
at all possible times Yrlthin Tg most first be considered. 



Stack 

Height 

S(t) 




37 



'\.8 llluatrated in Fig, 2.4, a group of planes arriring in the time 
interval (t, t ♦ dt), with t i T^., finds a stack of j( (g - a)t planss 
waiting to land. Now the landing rate is a R at any time t < T^, so 
that the delay for the grovcp of planes arriving during time dt is 

/y - ^ ^ ^ ~ ) 

< % ” 2.32 

o 6 6 “Tfv . 

Eq. 2,32 is only valid if the driay oT 1 - t, for if T > - t, all 

of the planes will not have landed by and we have to consider the 
change in landing rate at T^, Thus another limiting condition on t in 
Eq, 2,32, in addition to the ccoditlon 

O < ± j 

is that 

(£^)t i T.-t^ 

car 

o £ t i ^ X, , 

2.33 

The real limiting ccndltion on t in Eq, 2,32 therefore depends upon 
whether 

■''o 2 .% 



or 
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For the condition 



“T > t > T 



n. 



o 



( c«se A ) 
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the planes arriring in the interval dt will not have landed by T^. 
So the nxunber of planes still due to land before them at is 



Now it should bo clear that the maximum delay vriU occur for some 
group of planes arriving at or Just before Tj., where the stack takes 
on its maximum value and the landing rate remains small. The difficul- 
ties of finding the maximum near the end point of a linear equation 
of the form of So. 2.3^ need not be elaborated on. The method used 
here is empirical and the valv» of t obtained to maxlmlEe Er, 2.3S is 



1 5 cq - t ) oc g . 
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Because these planes land at a rate , 



O'; (T^-t ) -1- - t ^ - (t; - t ) c< ^ 

A 

- To ( i-cx) ^ -t ( » - e) . 



2,38 



t= . 



Thus 










Foft ^ ^ I ; Casff A * 
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2.40 



For the conditioia 

-t i T„ f To e) 

G ! 

the delay for any plane in the group is given slnply Eo, 2.32, and 
for T to bo a maximu m set 
t -T^. 

Thxis 

T - T f -L 

MAX - n 
for T o > ^ y. I 
2.5 Additional delay of aircraft carrier T , 

cLC 

It should be clear that the extra time needed by the aircraft car- 
rier to land aircraft vdien a low stack exists is simply 

*^ac ** ^s ~ *^r’ 2.2f2 

where 

T - N/r. 
r 

Thus T„- is the additional time required for the carrier to steam into 
ao 

the wind at high speeds and should be made a ndnimom for the most ef- 
ficient utHiiation of the landing platform. 

Case A: '> TjSl 't \ 

^ Tn. 

From I3qs. 2.4 and 2.42, 

•• To (!-<.) -~ r „ 0 - 1 ) . 2.W 



2.41 



} Case B. 
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Case B: T~ O- > ^ > | 

“Tn 

From EqSu 2o7 and 2.42, 




2.6 Summary of results. 

In the same manner, for T^, < the calculations split up into 
two portions, depending upon the time at which the stack finally clears. 
See FigSv 2,5 and 2,6, The ccunplete results are sxuimarized in Table 2, 
which gives the length of the stack period, the maximum stack size, the 
average delay, the maximum delay, and the aircraft carrier delay. 

From Table 2 the aircraft carrier delay T^^ is seen to be in general 
a function of the variables T^, T^, a , and £ , It is convenient to 
use a normalizing factor in preparing graphs for visual representation, 

T^ was chosen as the normalizing factor for convenience because in prac- 
tice it is usually fixed, A family of curves, for a fixed value of a , 
showing T^^Tj. vs the variable y = T^/T^ (to specify Tq) for 
several values of € as a parameter, were prepared. See Figs, 2,7, 2,8, 
and 2,9, Other curves could be prepared in a similar manner. 

The general implications of the numerical results may best be ap- 
preciated in tbrms of examples. Consider a carrier conducting landing 
operations with the ti'affic density £ approximately 1, Now let us 
assume that during the landing operations the landing rate has been 
suddenly reduced from the normal rate in such a manner that when a = 0,7, 
Tq = 0.4 Tj.; a - 0,5, T^ - 0,2 T^.; and a - 0,3, T^ - 0,2 T^., In 
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other words during landing operations the landing rate has been reduced 
in several instances so that on the average, say 0o8 of the time, land- 
ings were being made in an inefficient manner (poor pilot technique, 
inejqperienced plane handlers, etCo)<» Now from Figs, 2o7, 2.8, and 2»9 
we find that the additional delay to the carrier to be approximately 
O.i+Tj., That is, the carrier must spend almost half the time normally 
taken to land aircraft steaming at high speed using extra fuel in a 
most inefficient manner. Taken over a period of time this could amount 
to a sizeable wastage, to say nothing of the danger that might be 
present in enemy waters. 

Similarly, other values of the parameters may be used to ccxnpute 
new delays, and in addition the behaviour of the stack may be inves- 
tigated, i.e., maximum stack, average delays, etc. These nay be impor- 
tant in determining fuel reserves of the aircraft. 

The hazards in working with the traffic density €. approximately 
1 have been stressed throughout the paper. However, carrier landing 
operations are a unique application of queueing theory and it should be 
clear that the most efficient utilization of the landing platform will 
only be accomplished with £ approximately 1. And this can be accomplish- 
ed only vdien the pilots and other personnel are exercising the most pre- 
cise skill, Judgment, etc, during the landing operations. This is of 
course merely a statement of what common sense would dictate. 

2.7 Conclusions, 

For a simplified situation occuring in the landing of carrier air- 
craft, an investigation has been made into the variation of the average 
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aiin:raft waiting time, the time wasted by the carrier during landing 
operations, and the distribution of the number of aircraft in the stack 
for different initial conditions. The simplifications and assumptions 
necessary to ax^ive at a single solution make the problem somewhat arti- 
ficial, However, it is not completely unreasonable, because in actual 
practice the landing process may actually behave in some instances in 
a manner sljnilar to the process assumed. 

The practical verification and application of the results des- 
cribed in this paper depend upon the specification of the landing time 
distributions, the traffic density g , and the minimum safe landing 
interval. 

In comparing the present results with actual caiTier operations, 
many observations of this process must be made during actual landings 
so that the validity of the assumptions made may be verified. In prac- 
tice, then, a great deal of information is required to check the 
numerical model which has been used as a basis for the result pre- 
sented here. If these assumptions are not valid as based on observation- 
al data, the model discussed in the next section may be applicable. 
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2.8 Statistical IiapLlcatlons . 



The difficulties encoxintered In any analytical sdutlOT of a queue- 
ing process vdth regular arrivals and a respectatle distribution of ser- 
vice tlsBs Is apparent frcan Llndley's analysis (See 1.6). In particular, 
if the process is not in statistical equilibrium, the Individual waiting 
time distributions are difficult to obtain (Eq. 1.19). 

An alternative process is a quasl-emperlcal investigation using an 
electronic calculator as in the LOT report, or if that method is not 
available, the use of random numbers may be substlttcted . This is a well 
known standard alternative to a purely theoretical analysis, and its 
basic principle is as foUows. If we read off successive groups of, say, 
four or five numbers from tables of random nximbers, then these can be 
regarded as numbers, with four or five decimal places, drawn at random 
from a rectangular frer-uency distribution with limits 0 and 1. Now all 
continuous distribxition functions are themselves rectangularly distribu- 
ted between 0 and 1, so that the abscissa of any given distribution 
function corresponding to a value of the function taken fran a table of 
random numbers can bo regarded as a ouantity drawn at random, with the 
appropriate frequency, from the frequency distribution itself. Thus for 
a given distribution of landing time, for example a gamma distribution 
with given parameters, n and X say, we can use random nun±>ers to run 
off a series of independently chosen landing times. This series can then 
bo taken as the actual landing times of successive aircraft for a hypothe- 
tical aircraft carrier. Data for a largo number of hypothetical landing 
patterns can be obtained in this way, and the enqplrical distributions. 
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means, standard deviations, etc., for rariovis altematiro landing sys- 
tems can be worked out. 

Tilth this method of procedure it is clear that we have to decide 
on the basic frequency distribution of landing times beforehand, thou^ 
we could consider a limited range of types of distribution (rectangular, 
parabolic, etc . ) . 

In order to obtain a basic frequency distribution, observations Bust 
be made of actiial landings with the landing times recorded. The usual 
curve fitting techninues can then be employed. Once the proper distri- 
bution has been decided upon (this may not be easy) the method of maxi- 
mum likelihood may be used to find the unkncnm parameters from the 
sample data (Mood 10 ) . 

For conqputational ease, grouping of the frequency distribution could 
be utilised. From a practical point of view, a truncated distribiitlon 
may be more lealistic since the probability of extreme values are ef- 
fectly sero. .After the preliminary planning, the corresponding dis- 
tribution function for each type of curve decided upon must be calcula- 
ted to four decimal places. Four figure numbers are then read from the 
table of random numbers, and the value of x (landing tims) corresponding 
to the mid-points of the intervals into which the numbers fell are then 
recorded, . In this way series of landing times can be obtained readily, 
simply by taking successive groups of four figure numbers from the 
tables of randan numbers and writing down the corresponding value of 
X read from the calculated table of the appropriate distribution function. 

Now for any series of landing times the waiting time of each air- 
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craft and the carrier’s idle or waated tine at any time laay be calcu- 
lated by the usual set of recursive equations (see Fig. 1.1): 

'Y t nr + j? - n. 

If negative, that is simply the carrier’s idle time. Vdth a 

constant arrival interval, the computations are simple but tediovis. 

Finding the distribution of the actual stack is a little more com- 
plicated. The easiest method of computation seems to be to consider 
for how long in the Intervals (assiuned constant) between successive 
arrivals there are 0, 1, 2 aircraft waiting. For example, con- 

sider a constant departure interval of 1 unit. Consider five successive 
aircraft who have waits of 0, 0,5, 1.5, 2, and 1 units. Assume there 
is no stack when the first aircraft enters the landing pattern. Thus 
the first aircraft doesn’t wait, so there are no aircraft in the stack 
and it is empty for 1 unit, say. The second aircraft waits 0.5 units 
of the second interval, thus the stack is enqjty for the remaining 0,5 
tmits. The third aircraft >valts all of the third Interval and 0.5 units 
of the fourth interval j while the foxtrth aircraft waits aU of the fourth 
and fifth interval; and the fifth aircraft waits all of the fifth inter- 
val, Thus it follows that in the second interval there is one aircraft 
in the stack for 1 unit. In the fourth Interval there are two then caie 
aircraft waiting fer 0,5 units. It is not too difficult to carry out 
the arithmetical operations by using tables for visual representation 
as follows: 
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Waiting time 


Array 


Number of 


a/c in stack 
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0 1 


2 X 


0 




1 




0.5 




0.5 0.5 




1.5 




1 


0.5 


2.0 




0.5 


1 


1.0 









Th® waiting tines are grouped according to the arrival interval. The 
times during Tdiich 0, 1, , . . . , n aircraft are waiting during each inter- 
▼al are shown in the zdgbt hand table, and are obtained from successive 
diagonals of the center array. The figure appeaidng at the N. E. point 
of any diagonal shows how long the stack is occupied by the greatest 
number of aircraft in the corresponding interval. The latter nun4>er is 
given by the number of quantities in the diagonal. The difference be- 
tween this figure and the one immediately to its S« W. gives the length 
of time that one fewer aircraft are waiting, etc. Subtracting the figure 
appearing at the S. point of any diagonal from 1«0 gives the length 
of time that the stack is unoccupied. 

The above discussion explains how to obtain series of aircraft wait- 
ing times, carrier Idle time, and the stack sise. For any particular 
hypothetical carrier landing process one can calc\iLate the mean aircraft 
waiting timej the total carrier idle time; and the time distribution of 
the stack. By considering a series of such processes it is then possible 
to estimate what can be expected to happen in the long run, provided the 
observational landing data are vatlid. 
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